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Inclusive production of J/tj) mesons, especially at forward rapidities, is an important probe of 
small-® gluons in protons and nuclei. In this paper we re-evaluate the production cross sections in 
the Color Glass Condensate framework, where the process is described by a large x gluon from the 
probe splitting into a quark pair and eikonally interacting with the target proton or nucleus. Using 
a standard collinear gluon distribution for the probe and an up to date dipole cross section fitted 
to HERA data to describe the target we achieve a rather good description of the cross section in 
proton-proton collisions, although with a rather large normalization uncertainty. More importantly, 
we show that generalizing the dipole cross section to nuclei in the Glauber approach results in a 
nuclear suppression of J/if production that is much closer to the experimental data than claimed 
in previous literature. 

PACS numbers: 13.85.Ni, 14.40.Pq, 24.85.+p, 25.75.Cj 


I. INTRODUCTION 

Heavy quarks occupy an important place in QCD phe¬ 
nomenology. The mass of the charm quark is large 
enough that one should be able to describe its inter¬ 
actions using weak coupling. On the other hand, it is 
light enough to be a dynamic probe that is sensitive to 
momentum scales that characterize high energy nuclear 
collisions, such as the deconfinement temperature or the 
gluon saturation scale. Because of the clean experimen¬ 
tal signatures, vector meson production is a particularly 
important probe of heavy quark interactions. The “melt¬ 
ing” of the J/ip meson has traditionally been considered 
as one of the most important signals of the deconfinement 
transition in ultrarelativistic heavy ion collisions. Study¬ 
ing J/ip production in proton-proton and proton-nucleus 
collisions at high energies can also provide valuable in¬ 
sight into the physics of gluon saturation and strong color 
fields at small x, as well as into the energy loss of a high 
energy probe in nuclear matter. 

We will in this paper use the Color Glass Condensate 
(CGC) framework together with a simple color evapo¬ 
ration model to calculate cross sections for J/ip and T 
production at forward rapidities in proton-proton and 
proton-nucleus collisions at LHC energies. At forward 
rapidity the probe proton can be treated as a “dilute” 
collection of partons, and the CGC calculation proceeds 
according to a rather straightforward picture. The probe 
proton can be described by conventional collinear parton 
distributions in the so-called “hybrid” formalism, or by a 
corresponding “fc-r-factorized” approximation (which in 
fact strictly speaking violates fc^-factorization [ 1 ] since it 
involves several unintegrated gluon distributions). The 
physical picture in both cases is the same: a gluon from 
the probe can split into a quark-antiquark pair. As this 
system propagates through the target, it can pick up an 
eikonal phase described by an SU(3) color matrix, the 
“Wilson line” from the target. This leads to an expres¬ 
sion for the cc pair production cross section (in either 
singlet or octet color states) in terms of Wilson line cor¬ 


relators describing the target. The same Wilson line cor¬ 
relators, the simplest one of them being the “dipole cross 
section”, appear in calculations of e.g. total DIS cross 
sections [2-4], single [5-8] and double inclusive [9-12] 
particle production in proton-proton and proton-nucleus 
collisions, diffractive DIS [13, 14] and initial state for the 
liydrodynamical modeling of a heavy ion collision [15-17] , 
making the framework very broadly applicable to many 
different processes. 

Several LHC experiments have performed detailed 
measurements of inclusive J/ip production in a broad 
range in rapidity, particularly in proton-proton [18-22] 
but also in proton-nucleus [23, 24] collisions. The results 
on the “cold nuclear matter suppression”, i.e. the sup¬ 
pression of the cross section in proton-nucleus relative 
to proton-proton collisions, have then been compared to 
several theoretical predictions. In particular, the nuclear 
suppression predicted by the CGC calculation in [25] has 
been much stronger than observed in the data. The main 
purpose of this paper is to explore to what extent this 
disagreement stems from the framework itself, and how 
much is due to approximations concerning nuclear geom¬ 
etry that are necessary for relating the CGC description 
of the proton to that of the nucleus. Indeed, in the case 
of single inclusive light hadron production [8], the dis¬ 
agreement with early CGC nuclear suppression calcula¬ 
tions [26] and LHC data turned out to be mostly due 
to the nuclear geometry effects. We will in this paper 
calculate cross sections for inclusive (prompt) forward 
J/ip production in proton-proton and proton-nucleus col¬ 
lisions using the collinear “hybrid” framework also used 
in Ref. [25]. 

We use here the dipole cross sections obtained in 
Ref. [8] by performing a fit to HERA data using the run¬ 
ning coupling BK equation. These have several advan¬ 
tages over the ones used in previous works. Firstly the 
dipole cross section in a nucleus is related to the one in 
a proton taking into account the fluctuating positions of 
the nucleons in the nucleus via a Glauber model, without 
any additional parameters beside the standard Woods- 
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Saxon nuclear density. This leads, unlike many previous 
CGC calculations, to nuclear modification ratios that ex¬ 
plicitly approach one in the high transverse momentum 
limit, at all collision energies. Physically this should be 
expected from the validity of collinear factorized pQCD 
in this regime. Secondly, the dipole cross section from 
Ref. [8] corresponds to an explicitly positive unintegrated 
gluon distribution when Fourier-transformed to momen¬ 
tum space. Thirdly, our calculation, like that of [8], takes 
consistently into account also the proton transverse area 
from the fit to HERA data. This completely fixes the 
normalization of the LO calculation at the quark level 
without any additional parameters to describe the pro¬ 
ton size 1 . The use of also the proton transverse area 
from the DIS fits and the way of extending the dipole 
cross section from a proton to a nucleus are quantita¬ 
tively the most important differences between our work 
and that of [25]. 

This paper is organized as follows. We will first briefly 
describe the formalism used in the calculation in Sec. II 
and the used dipole correlators in Sec. III. We will then 
present our results for cross sections in proton-proton 
and proton-nucleus collisions in Sec. IV and finish with 
a discussion and outlook for future work in this direction 
in Sec. V. 


II. FORMALISM 


We will here use the simple color evaporation model 
(CEM), where a fixed fraction of all cc pairs produced be¬ 
low the D-meson threshold are assumed to become J/tj) 
mesons. The quark pair can be produced in either the 
singlet or octet representation, the additional color is as¬ 
sumed to be “evaporated” away in the form of a soft 
gluon. In the CEM model, the differential cross section 
for J /■0 production therefore reads 


_ p [ D i n ,r2 dg- ee 

d 2 Pj_ dY J/ ^ J im 2 d 2 Pj_ dY dM 2 ’ 


(1) 


where P± and Y are the transverse momentum and ra¬ 
pidity of the produced J/ip, d 2 P±dYdM 2 1S the cross sec¬ 
tion for the production of a cc pair with transverse mo¬ 
mentum -P_l, invariant mass M and rapidity V, m c is the 
charm quark mass and mjj = 1.864 GeV is the D meson 
mass. Here Fj/i/, is a nonperturbative quantity repre¬ 
senting the probability that a cc pair has to form a 
In the literature [27] this constant is taken to have dif¬ 
ferent values depending on the c-quark mass, the parton 
distribution and the factorization scale. We will here use 


1 Note, however, that the normalization is expected to be quite 
different at NLO, and that for J/ip in the color evaporation model 
the normalization is given by an additional phenomenological 
parameter. 


a constant value Fj/ip = 0.06. This is somewhat larger 
than the typical values in perturbative NLO calculations, 
but it should be kept in mind that our calculation is at 
the LO level and is therefore likely to underestimate the 
normalization of the cc cross section. The same model 
can be used to compute the cross section for T produc¬ 
tion by simply replacing m c with the bottom quark mass, 
Mr) with the B meson mass Mb = 5.280 GeV and Fj/^ 
with Ft that we take equal to 0.04 here. One should 
note that since the 6-quark mass is already significantly 
higher than the saturation scale, bottom quark produc¬ 
tion is not in the natural regime of validity of the CGC 
framework. We will in this paper however calculate also 
T cross sections for comparison. 

The formalism for analyzing gluon and quark pair pro¬ 
duction in the dilute-dense limit of the CGC formalism 
has been worked out in great detail already some time 
ago [28, 29] (see also [30]) and applied in many calcula¬ 
tions such as [25, 31-33]. It can be obtained from the 
observation that an incoming gluon from the projectile 
can split into a quark-antiquark pair either before or af¬ 
ter the interaction with the target. These partons are 
then assumed to eikonally interact with the target, pick¬ 
ing up a Wilson line factor in either the adjoint or the 
fundamental representation, depending on the particle. 
We will here use the collinear approximation, where the 
incoming gluon is assumed to have zero transverse mo¬ 
mentum. Calculated in this way the cross section for 
the production of a cc pair with transverse momenta p Y 
and q Y and rapidities y p and y q reads, in the large--/V c 
limit [25], 


_ da C c _ _ a 2 N c 1 l "B con (p ± +q ± ,k±) 

d 2 p ± d 2 q ± dy p dy q 8t r 2 d A (2? t) 2 J (p ± + q ± ) 2 

k± 

X in jJP-L +<7±> fe J-) x lGp(Xi,Q 2 ) , (2) 

x 2 


where J k ^ = J d 2 k±/(2n) 2 , and d A = N 2 — 1 is the 
dimension of the adjoint representation of SU(lVc). The 
longitudinal momentum fractions of the gluons coming 
from the projectile and the target, X\ and x 2 , are given 

by 


Xl,2 


VP± 2 + M 2 ±Y 


( 3 ) 


The “hard matrix element” c, co n is the sum of 3 terms, 


77 _ 7799.99 I ^99,9 I 7:9,9 

1 'coll — 1 'coll ' “coll “ r “coll ! 


( 4 ) 


corresponding to a quark pair scattering off the target 
( qq 1 qq ), the gluon scattering before splitting into a quark 
pair ( g , g ) and an interference between the two. The 
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explicit expressions are 
8 p + q + 


■=" 39.99 

^coll 


^qq,g 

^coll 


^9,9 

‘“'coll 


(p + + q + ) 2 (a ± 2 + ml ) 2 

m 2 , b + ) 2 + (? + ) 2 a 2 

TTl r - - -;---- a L 


(p+ +q+) 2 
16 


(p + q) 2 {a± 2 + ml) 

2 , (p + ) 2 + (q + ) 2 


(p+ + q+) 3 


a± ■ (p + q_ l - q + p±) 


(jp + q ) 4 


(.p + q ) 2 - 


(p+ + q+) 2 


(p + q±- q + p± f 


(5) 



with a l = — fcj_. 

The propagation of a quark-antiquark pair through the 
color field of the target is described by the function 


FIG. 1: Differential J/ij) cross section as a function of Y in 
pp collisions. Data from Refs. [19, 22]. 


/ AT 1 2 

d 2 h±-^ S Y (k ± ) S Y (l x - k ± ) , 

( 6 ) 

where b± is the impact parameter. All the information 
about the target is contained in the function S Y (k±), 
which is the fundamental representation dipole correlator 
in the color field of the target: 

S Y (k ± ) = J d 2 r ± e ik ^S Y (r ± ) , (7) 

with 

SA*± ~ V.l) = ^-(TrU\ Xl _)U(y ± )) , (8) 



where U (x_l) is a fundamental representation Wilson line 
in the target color field. Our calculation is a leading order 
one, so for consistency we should use a leading order 
collinear gluon distribution to describe the probe. In this 
work we use the MSTW 2008 [34] LO parametrization for 
G p {xi,Q 2 ). 


III. DIPOLE CORRELATOR 


For the dipole correlator S Y (k±) we use the MV e - 
parametrization from Ref. [8] . This is obtained by fitting 
the combined inclusive HERA DIS cross section data [36] 
in the region Q 2 < 50 GeV 2 and x < 0.01 using a dipole 
cross section, whose ^-dependence is obtained with the 
running coupling BK equation [37-39]. The dipole am¬ 
plitude at the initial rapidity Xq = 0.01 is given by the 
parametrization 


S Y= in = exp 

XQ 


r± 2 Q, 


sO 


In 


( 1 

VIglIAqcd 



(9) 

The impact parameter profile in the proton is assumed to 
factorize from the dimensionless dipole amplitude, and is 


FIG. 2: Differential T cross section as a function of Y in pp 
collisions. Data from Refs. [22, 35]. 


parametrized by a constant 

J d 2 b± -A 2 ° , (10) 

so that the distribution needed for quark pair production 
is given by 

4>l% 9 (l±,k ± ) = ^ S Y (fcjJ S Y (l ± - fcx). (11) 

Using the expression 


( 33 - 2 A)log(^) 

for the running coupling constant, the fit gives a very 
good description y 2 /d.o.f = 1.15 of the data with the 
parameters Q 2 0 = 0.060 GeV 2 , C 2 = 7.2, e c = 18.9 and 
CTo/ 2 = 16.36 mb. One advantage of the parametriza¬ 
tion [8] over the AAMQS [2] one used in [25] is that 
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FIG. 3: Differential J/ip cross section as a function of P± in 
pp collisions (2.5 < Y < 4). Data from Refs. [19, 22]. 
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FIG. 4: Differential T cross section as a function of P± in pp 
collisions (2.5 < Y < 4). Data from Refs. [22, 35]. 


the Fourier transform of Sy’=i n J- ( r J-) is positive definite 

IQ 

even at the initial rapidity, enabling a more natural phys¬ 
ical interpretation as an unintegrated gluon distribution. 
We also used the MV and MV 7 parametrizations from 
Ref. [8] for the initial condition of the BK equation and 
found that the effect of such a change is small compared 
with the other uncertainties involved in the calculation. 

For the nucleus we use, again following [8], a BK- 
evolved dipole correlator with the initial condition 


with d = 0.54fm and Ra = (1.12A 1 / 3 — 0.86A -1 / 3 ) fm, 
is the standard Woods-Saxon transverse thickness func¬ 
tion of the nucleus, normalized to unity. To calculate 
cross sections with nuclear targets we explicitly integrate 
over the impact parameter, treating the dilute edge of 
the nucleus (where the saturation scale of the nucleus 
falls below the proton saturation scale) as in Ref. [8] . In 
practice this amounts to using the proton-proton result 
scaled such that the nuclear modification factor R p a , de¬ 
fined in Eq. (15), is unity in the dilute edge of the nucleus. 
This “Glauber model” form can be derived by assuming 
that at the initial rapidity, a high energy probe scatters 
independently off the nucleons in the nucleus. The av¬ 
erage over the fluctuating positions of the nucleons can 
then be performed analytically. At small r the integral of 
Eq. (13) over b explicitly approaches A times the dipole- 
proton cross section (9). Thus the form (13) explicitly 
leads to nuclear modification ratios R p a (see definition 
Eq. (15)) that, at all collision energies, approach unity 
at small distances, i.e. large momenta. We emphasize 
that besides the standard nuclear density Ta (b) there 
are no additional parameters in passing from the proton, 
Eq. (9), to the nucleus in Eq. (13). Also the transverse 
area of the proton cro/2 is consistently taken to be the 
same in proton DIS and in hadronic collisions, i.e. the 
identification (10) is also used in the pp —► J/ip + X cross 
section. 

In our calculations the nuclear modification only de¬ 
pends on parameters fit to proton data and on the stan¬ 
dard nuclear thickness function. In particular, no inde¬ 
pendent nuclear saturation scale with an unknown value 
is needed, in contrast to Ref. [25]. The practical result of 
this is that the effective typical nuclear saturation scale 
from Eq. (13) is smaller than the A 1//3 -scaling from the 
proton assumed in Ref. [25]. On the other hand, we treat 
separately the area occupied by the small -x gluons in the 
proton, <Jo/2, and the total inelastic cross section, <Ti ne i. 
Thus here the J/ip cross section in proton-proton colli¬ 
sions is proportional to (Jo/2 which is smaller than the 
proton area assumed in Ref. [25]. The net effect of the 
smaller nuclear Q s and the smaller yield in proton-proton 
collisions is that one obtains a nuclear suppression ra¬ 
tio R p a that is closer to one around P± ~ Q s than in 
Ref. [25], and simultaneously approaches unity at large 

P±- 


IV. RESULTS 


s li„i(u,&i) = exp 


-AT A (b±) 


(Jo r± 2 Q 


2n2 

s0 


2 4 


x In 


1 


|gl|Aqcd 


+ e c • e 


where 


T A (b±) = j dz- 


n 


\J b± 2 +z 2 -RA 


d 


(13) 


(14) 


In this section we will present the results of our calcu¬ 
lation and show a comparison with LHC data from the 
ALICE and LHCb experiments in pp and pA collisions 
at center of mass energies of 7 and 5 TeV respectively. 
It should be noted that both for pp and pA collisions, 
LHCb gives separately the cross section for prompt J /ip 
production and the contribution from b decays (which is 
not included in our calculation) while ALICE data is only 
for inclusive 3/ip production. Therefore, for a consistent 
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FIG. 5: Differential J/ip cross section as a function of Y in FIG. 6: Differential T cross section as a function of Y in pPb 
pPb collisions. Data from Refs. [23, 24]. collisions. Data from Ref. [40]. 


comparison, we subtract the b decay fraction measured 
by LHCb from the ALICE data to extract the prompt 
contribution. 

In the uncertainty band we include the variation of 
the charm quark mass between 1.2 and 1.5 GeV (in the 
case of J/ip production) or of the bottom quark mass 
between 4.5 and 4.8 GeV (in the case of T production) 
and of the factorization scale between 2M± and M±/2 
where M± = \JM 2 + P± 2 and M is the qq pair’s invari¬ 
ant mass. When we show quantities integrated over the 
transverse momentum of the produced meson, we inte¬ 
grate up to P±_ = 15 GeV. Our results do not depend sig¬ 
nificantly on this choice since the cross section is quickly 
decreasing at large P±. 
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FIG. 7: Differential J/ip cross section as a function of P± in 
pPb collisions (1.5 <Y< 4). Data from Ref. [24]. 

A. Proton-proton collisions 


We show our results for the Pj_-integrated cross section 
d< ^y V ' at a center of mass energy ^s = 7 TeV in Fig. 1, 
compared with data from ALICE [22] and LHCb [19]. 
We see that the shape of the data is quite well described, 
but the normalization uncertainty is quite large. Similar 
conclusions can be drawn for T production, as can be 
seen from Fig. 2 where we compare our results with AL¬ 
ICE [22] and LHCb [35] data, although we observe that 
the uncertainty band of our calculation is narrower than 
for J/ip production. 

In Fig. 3 we show the comparison of our results for the 
J/ip P± spectrum with ALICE and LHCb data. Again 
the agreement with the data is quite good except at very 
large transverse momenta where the CGC description is 
not expected to be very accurate. For T production, 
as shown in Fig. 4, the agreement is also quite good at 
higher P±, but this leading order framework does not 
describe the increase in the mean P± of the T compared 
to the J/ip. 


B. Proton-lead collisions 

In Fig. 5 we show the differential J/ip cross section as 
a function of the rapidity in pA collisions at a center of 
mass energy of 5 TeV and compare with ALICE [23] and 
LHCb [24] data. As in the pp case, the slope of the data is 
compatible with the one predicted by our calculation but 
the normalization uncertainty is large. In Fig. 6 we show 
the same quantity for T meson production and compare 
our results with data from the ALICE Collaboration [40] . 
Here we can see that the normalization uncertainty of our 
calculation is again smaller than for J/ip production, but 
there are only two data points to compare with, which 
limits the scope of this comparison. 

In Fig. 7 we show the differential cross section as a 
function of the transverse momentum of the J/ip in pA 
collisions at a center of mass energy of 5 TeV and com¬ 
pare with LHCb [24] data. The agreement with data 
is good but it is not possible to see if our calculation 
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RpPb (J/lp) 



FIG. 8: Nuclear modification factor for ,J/■i/i production as a 
function of Y. Data from Refs. [23, 24] 


RpPb (J/ip) 



FIG. 10: Nuclear modification factor for J/-i/> production as a 
function of P± (2 < Y < 3.5). Data from Ref. [42]. 


RpPb (T) 



RpAu (J/lp) 



FIG. 9: Nuclear modification factor for T production as a 
function of Y. Data from Refs. [40, 41]. 


overestimates the cross section at large P± like it does 
for pp collisions since there is only one experimental bin 
between 8.5 and 14 GeV. 

The normalization uncertainty from the charm quark 
mass and the factorization scale is quite large in the cross 
section itself. The modification of the production mecha¬ 
nism due to the stronger color fields in nuclei is, however, 
very little affected by these factors. Indeed, if the nuclear 
geometry is properly treated within the Glauber picture, 
the CGC framework provides a very robust prediction for 
the nuclear modification factor, defined as 


_ i da/d*P ± dY\ pA 
PA A da/d*P ± dY\ pp • 


(15) 


This ratio has been extracted from the proton-proton and 
proton-nucleus cross sections by both the ALICE and 
LHCb collaborations. In Fig. 8 we show it in the case of 
J/ij) production as a function of the rapidity and compare 
with experimental data. We observe that our results are 


FIG. 11: Nuclear modification factor for J/i)> production at 
RHIC energies as a function of P ± (1.2 <Y< 2.2). PHENIX 
data from Ref. [43] 


close to the data when taking into account all sources 
of uncertainty. In particular, as already discussed in the 
Introduction, our calculation is much closer to the data 
than early CGC calculations in Ref. [25]. The nuclear 
modification factor for T production was also measured 
by the ALICE [40] and LHCb [41] collaborations, in two 
and one rapidity bins respectively. We show the com¬ 
parison of these data with our results on Fig. 9. Here 
no firm conclusions can be drawn at the moment since 
experimental uncertainties are quite large. 

Finally we look at the the behavior of R p a as a func¬ 
tion of P±. Figure 10 shows the behavior of R p a as 
a function of P± for J/ip production, compared to the 
measurement by the ALICE collaboration [42]. Except 
at very low P ±, the agreement of our calculation with this 
measurement is quite good, and the normalization uncer¬ 
tainty has almost completely canceled out. In Fig. 11 we 
also compare our calculation to PHENIX results [43] at 
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lower energies. While our calculation displays a quite 
strong “Cronin” peak that is not seen in the data, also 
here our results are closer to data than the previous CGC 
calculations [25]. The nuclear modification factor results 
at LHC energies, Figures 8 and 10, are the main result 
of this paper. 

V. DISCUSSION AND OUTLOOK 

At high enough transverse momentum, particle pro¬ 
duction in QCD processes should be well described by 
collinear factorization. Indeed many comparisons of NLO 
calculations using nuclear parton distributions have been 
performed in the literature, see e.g. Ref. [44], In gen¬ 
eral it seems that the shadowing present in the standard 
EPS09 [45] nuclear gluon distribution tends to underes¬ 
timate the nuclear suppression observed experimentally 
in proton-nucleus collisions at lower transverse momenta. 
The collinear NLO calculation is in fact not that far from 
the physical picture of our work here. A leading order 
collinear calculation of heavy quark pair production re¬ 
sults in a delta-function P^-distribution for the pair. To 
get a physically meaningful result for the differential cross 
section, one must go to NLO, where an additional radi¬ 
ated gluon can carry the recoil transverse momentum. 
In a standard NLO calculation, one integrates over the 
phase space of this gluon. At high energy this additional 
phase space integration can contribute a large logarithm 
~ Ins. The BK equation is, in this context, a resum¬ 
mation of these logarithmically enhanced terms to all or¬ 
ders in a s lns. The large nuclear suppression in the data 
could be taken as an indication that the nonlinear effects 
included in BK indeed are significant for the J/ip. 

In “cold nuclear matter” (CNM) energy loss calcula¬ 
tions such as Ref. [46], the physical argument behind the 
calculation can also be described in a language that is 
very similar to the one used here. The picture in the tar¬ 
get rest frame is that of an incoming gluon splitting into 
a quark-antiquark pair. The pair can then lose energy as 
it propagates through the target nucleus, leading to a nu¬ 
clear suppression of the cross section. This picture does 
not rely on the eikonal approximation, making it possible 
to extend the formalism to higher xa and transverse mo¬ 
menta. The practical implementation of the CNM model 
relies on a factorization of the quark pair production cross 
section and subsequent propagation through the nucleus. 
This factorization becomes questionable at very high en¬ 
ergies in the target rest frame (very high s and forward 
rapidities), when the typical target transverse momen¬ 
tum scale Q s is of the order of the charm quark mass. In 
this limit the timescales of the quark pair formation and 


of the gluon radiation leading to the energy loss are not 
well separated. Indeed, in the CGC calculation the two 
are not factorized, but the pair production and propaga¬ 
tion through the target nucleus are described coherently 
in the same formalism. 

There is a normalization uncertainty in our calculation 
from the fact that we are using the dipole cross section [8] 
that was obtained for light flavors only. Requiring sepa¬ 
rately a good description of the charm DIS cross section 
Fg would require additional parameters in the fit, e.g. 
as a separate cr 0 for heavy quarks as in the AAMQS [2] 
parametrization. This would improve our confidence in 
the normalization of the charm production cross section, 
but would be a significant work outside the scope of this 
paper. 

We also used the very crude color evaporation 
hadronization model to obtain the J/ip cross section from 
the quark pairs. Our estimate for the normalization un¬ 
certainty was rather conservative, as we varied the quark 
mass without adjusting the value of the color evaporation 
model constant Fj together with it. The description of 
the transition from a quark pair to a J/ip bound state has 
seen many recent theoretical advances in terms of non- 
relativistic QCD (NRQCD) [47, 48]. This has recently 
been successfidly combined [49, 50] with the small-a; for¬ 
mulation for quark pair production. It would be inter¬ 
esting to combine the newer proton and nucleus dipole 
cross sections used here with the NRQCD hadronization 
picture. This is, however, left for future work. 

Both of the uncertainties discussed above affect more 
seriously the separate proton-proton and proton-nucleus 
cross sections. They should mostly cancel in the nuclear 
modification factor R p a- We leave a more detailed inves¬ 
tigation of the separate cross sections to future work. The 
main purpose of this paper has been to focus on the im¬ 
portance of a proper treatment of the nuclear geometry, 
including the use of a proton transverse area consistent 
with the DIS fits used to obtain the dipole cross section. 
We have shown that when this is done using the dipole 
cross section parametrization obtained in Ref. [8] , the re¬ 
sulting nuclear modification ratio is much closer to the 
experimental data than in previous CGC estimates. 
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